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Abstract. Mean motion resonances (MMRs) can lead either to chaotic or regular motion. We report on a numerical experiment 
showing that even in one of the most chaotic regions of the Solar System - the region of the giant planets, there are numerous 
bands where MMRs can stabilize orbits of small bodies in a time span comparable to their lifetimes. Two types of temporary 
stabilization were observed: short period (~ 10"* years) when a body was in a MMR with only one planet and long period (over 
10^ years) when a body is located in overlapping MMRs with two or three planets. 

The experiment showed that the Main Belt region can be enriched by cometary material in its pre-active state due to temporary 
resonant interactions between small bodies and giant planets. 
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1. Introduction 

Commensurabilities in a mean motion are very common in the 
Solar System. Resonant relations in the Main Belt have been 
known for nearly one and a half century - gaps in the dis- 
tribution of asteroids were first linked to mean motion reso- 
nances (MMRs) by Kirkwood (1867). The influence of com- 
mensurabilities on the dynamics of small celestial bodies can 
be found in every region of the Solar System: from the Kuiper 



Belt (Nesvorny & Roig EOOOl l to the vicinity of the Venus orbit 
(Milani & Baccili 1999i 

Resonant relations of small bodies in the region among 
giant planets have been studied for the last three decades. 
Long time orbital stability was found close to the L4 and L5 
Lagrangian points of Jupiter (Everhart 1973), in the narrow 
bands centered on 7.02 and 7.54 AU (FrankHn et al. 1989 1 and 
in the region between 24 - 27 AU (Holmanjl997i. Recent pa- 



pers describe the local (Grazier et al. ll999> or global (Robutel 
& Laskar l2001t dynamics in the region of the giant planets. 

In the cited works the exploration of small body dynam- 
ics was used to find stable regions of the Solar System where 
particles can survive on Gyr timescales on circular or quasi- 
circular orbits. However, the aim of this work is to investigate 
the dynamics of real particles in order to find the possible paths 
of their evolution and to learn about the role of MMRs as the 
mechanism of the orbital temporary stabilization process. 
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2. Method of computations 

In most of the previous papers on the subject, the results were 
based on the integrations of the equations of motion for fic- 
titious test particles, with initial orbital elements chosen 'ad- 
hoc'. In this paper, on the contrary, we simulate the evolu- 
tion of the orbital elements of real Centaurs and their clones. 
We chose orbits of all objects situated entirely inside the gi- 
ant planets region: 1998 SG35, 1999 UG5, 2000 fiCgg, 2000 
GM137 and 2001 BL41. All these objects are unnumbered, of- 
ten single-oppositional which means that their elements stored 
in public databases can be approximated. That is why we de- 
cided to re-calculate their orbits using all available observa- 
tions. Unfortunately, the data did not allow us to acquire the 
elements of 2000 GMuj with a reasonable precision, so we 
omitted this object in our computations. 

To compute the orbital elements of these real Centaurs at 
the initial epoch of our simulations, we took the observations 
from the files available at Minor Planet Center in Cambridge, 
USA. The sets of apparitions covered a very short time in- 
terval of about a month to a month and a half for the 1998 
5G35, 2000 ECgs and 2001 BL41, and over 4 months for the 
1999 UG5. Observations were selected according to the math- 
ematically objective criteria elaborated by Bielicki and Sitarski 
( fl991j . assuming a normal distribution of their random errors. 
Only 1 observation of the 1998 SG^s and 4 observations of 
the 1999 UG5 were rejected, the rest were used in the orbit im- 
provement process. Calculations were executed using Sitarski's 
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Table 1. Orbital elements and their mean errors of the four bodies derived from the observational data on the epoch 1998 07 06 
ET 



Object 


q±dq [AU] 


e±de 


i±di [deg] 


RMS [arcsec] (number of obs.) 


tL)±da> [deg] 


n±dQ. [deg] 


T±dT [ET] 


1998 5 G35 


5.84069270+0.00028536 


0.30587190+0.00001753 


15.62341+0.00011 


0.38 (49 observations) 


337.90443+0.00385 


172.49261+0.00008 


2008 02 10.31886 + 0.9552 


1999 UGs 


7.46349029+0.00027588 


0.41436520+0.00006159 


5.58734+0.00006 


0.61 (165 observations) 


288.36650+0.01335 


87.02353+0.00058 


1998 11 29.32162 + 0.21010 


2000 ECgg 


5.79188701+0.00393862 


0.45879307+0.00049997 


4.35728+0.00068 


0.48 (38 observations) 


163.21158+0.19737 


172.62330+0.00091 


2015 05 05.81412 + 15.12784 


2001 BL41 


6.84577522+0.21428445 


0.30172882+0.02988279 


12.50754+0.16842 


0.49 (29 observations) 


131.62288+6.14676 


280.58723+0.18671 


1998 01 20.61111 + 68.56299 



orbital programme package (Sitarski fDTTIIDT^FTW?! . Table 
1 gives the elements acquired for the bodies. 

2.1. Particle cloning 

These small bodies move in the giant planets region, with the 
perihelion distance close to Jupiter orbit and with the aphelion 
between the orbits of Saturn and Uranus. Orbits of the par- 
ticles are chaotic due to strong planetary perturbations, their 
Lyapunov time is short and lasts about 1000 years. 

As an example, let us take two orbits of bodies moving in 
the gravitational field of the 4 planets - their 5 orbital elements 
are identical and they only differ in the starting semi-major axes 
by 1.5 km (10"** AU). Fig. lA shows the eccentricity variations 
as a function of semi-major axis for two orbits - the nominal 
(crosses) and the changed one (circles) in a time span of only 
1 1,000 years. Fig. IB shows the logarithm of the divergence of 
semi-major axes between these orbits. The divergence is close 
to 10"^ AU in the first 3,000 years of the evolution but after 
this time it changes rapidly and both orbits start to evolve dif- 
ferently (see the huge leap of values of the logarithm on Fig. 
IB). This means that even very small changes of one of the 
elements make the particle's dynamics quite different after a 
small time span. To avoid this problem, an ensemble of possi- 
ble orbital elements is often used in such cases - the object's 
dynamical evolution is described in terms of probability. 

Usually sets of elements are acquired by varying initial or- 
bital elements within a reasonably small range, in most cases 
in the range of their mean errors. This procedure forms differ- 
ent orbital element sets but they cannot be treated as element 
sets of one celestial body. In fact they represent similar orbits 
of different bodies since these orbits do not fit the observations 
well. The sets of elements, used in the modelling and described 
in this paper, were generated in a different way. We used the 
method (found by Sitarski 1998 1 for creation of orbits which 
represent the observations well. This method makes it possi- 
ble to produce any number of orbital element sets and all of 
them fit the observations. Hence, these orbits can be treated as 
possible orbits of one celestial body. 

The dynamics of small bodies was modelled with the use 
of 404 orbital element sets created by the method mentioned 
above: for every real body 100 additional orbits were created. 
Table 2 shows the boundary values of all six elements for each 
real object. 




6 8 10 12 14 16 18 20 

Semi-major Axis a [AU] 




2000 4000 6000 8000 10000 
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Fig. 1. A. The phase space of the eccentricity variations in a 
function of semi-major axis for two nearly identical orbits: the 
reference (crosses) and the changed one (circles). They differ 
only in the starting semi-major axes by a factor of 10"^. The 
integration time is 1.1 ■ 10'* years. B. Time variations of the 
logarithm of the divergence in semi-major axes between orbits 
presented in the A. The leap near T=3000 is caused by two 
close encounters with Saturn of both objects. 

3. Numerical Model 

Dynamical evolution was investigated by integrating a 6-body 
problem, with the central mass, 4 giant planets from Jupiter 
to Neptune, and a massless particle. We analyzed 404 sets of 
initial conditions, each of them consisting of positions and 
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Table 2. Dispersion of objects' orbital elements - boundary values chosen from the set of 100 additional orbits, separately for 
each single element 



Object 


T 


q 


e 


w 


Q 


i 


1998 5 G35 


2008 02 10.55929 


5.84125774 


0.30592525 


337.91370 


172.49281 


15.62380 




2008 02 10.07866 


5.83989244 


0.30583258 


337.89299 


172.49240 


15.62309 


1999 UG5 


1998 11 29.85783 


7.46422563 


0.41450678 


288.40122 


87.02481 


5.58749 




1998 11 28.91046 


7.46296851 


0.41424401 


288.34168 


87.02197 


5.58720 


2000 ECgs 


2015 06 12.06479 


5.80080946 


0.46019067 


163.70350 


172.62571 


4.35859 




2015 03 24.73689 


5.78133382 


0.45761563 


162.65756 


172.62160 


4.35549 


2001 SL41 


1998 11 18.53404 


7.34292453 


0.40165728 


152.56526 


281.04189 


12.98724 




1997 10 26.14988 


6.23554911 


0.24252707 


119.25777 


280.07133 


12.10321 



velocities for the planets (the same values in all cases) and 
for the massless particle. The orbital elements for the plan- 
ets were taken from the Bretagnon's planetary theory. For the 
massless particle, we used the 4 orbits of real Centaurs and 
400 orbits cloned using Sitarski's method (100 clones for each 
real body). The initial positions and velocities were reduced to 
the barycentre of the inner Solar System. The sets were inte- 
grated separately 200,000 years forwards using the 15-th order 
RADAU integrator (Everhart [D7?l I1985> with an eiTor toler- 
ance set to 10"'^. This integrator adjusts the step size to main- 
tain the accuracy for all objects taking part in the process. 

Locations of mean motion resonances were simply calcu- 
lated with the use of Kepler's third law. The critical angle of 
the resonance was defined as follows: 

cr — ki ■ Ap — k2 ■ A + (k2 — ki) ■ (I) 
where ki , k2 are integers, Ap and A are the mean longitudes of a 
planet and a test particle respectively, and a> denotes the longi- 
tude of the perihelion of test orbit. During the experiment only 
two-body MMRs were examined. 

We accepted a special method for a resonance description. 
In our notation, a resonance 2/3 means that a planet takes 2 
revolutions around the barycentre in the time in which a small 
body takes 3. 

4. The influence of MMRs on small body dynamics 

The experiment showed that MMRs can temporarily stabilize 
numerous elliptic orbits in the Jupiter - Uranus region in peri- 
ods of time comparable to their typical lifetimes (i.e., the time 
over which non-resonant particles can survive in that region). 

Two types of resonant behaviour were observed. The first 
one - when a small body is in a two-body MMR with only 
one planet. This relation usually lasts a short time, from several 
hundred up to 35,000 years. There were 120 particles of this 
resonant type. We only noticed two particles in a resonance 
with one planet for a longer time span. 

One of these particles, corresponding to one of the 1998 
SG35 clones, is shown in Fig. 2. After 1 18,000 years of chaotic 
evolution the particle enters the 2/3 resonance with Jupiter and 
remains in this relation to the end of the integration process. 
Fig. 2A presents variations of the semi-major axis as a func- 
tion of time. The critical angle usually librates about 0° when a 
body is located in this particular resonance (see Fig. 2B). Even 
high and growing amplitude (over 150°) of the librations is in- 



capable of destabilizing the dynamics of the particles in a time 
span of over 80,000 years. 

Similarly to this Hilda type orbit, nearly all other orbits 
remaining in resonance with only one planet behave like the 
Toro class orbits (see Fig. 3). This class defined by Milani 
J 19891 1 consists of Earth-crossing objects which avoid close 
approaches by entering a MMR. Most particles remaining in 
the resonant relation with Jupiter or Saturn cross their orbits 
in ecliptic projections but cannot approach them. The time of 
such resonant relations is smaller than for Toro like orbits (with 
some exceptions) but this can be caused by stronger perturba- 
tions acting on a body in the giant planet region. 

The second type of resonant behaviour is characterized by 
much longer resonant interactions up to 160,000 years in our 
sample. Test particles are in two-body MMRs with 2 or 3 plan- 
ets (see Fig. 4). These overlapping resonant relations can sta- 
bilize orbits far more effectively - for over 10^ years. Fig. 4 
presents the evolution of the MMRs critical angles for one 
of these orbits. When two MMRs overlap each other, the be- 
haviour of the critical angle of the dominant resonance starts 
to vary (see Fig. 4C) - the critical angle begins to circulate 
with an overlapped libration of amplitude almost constant. 
Irregularities of the rotation and of the amplitude of libration 
do not mean that the body will leave the MMRs in the near fu- 
ture but the possibility that a small body will approach a planet 
increases with time. Only 9 to 15 % of orbits were stabilized on 
time spans longer than 35,000 years (according to the object). 

Fig. 4 shows that overlapping resonances can also protect 
the body from close approaches to planets. The disturbing res- 
onance does not necessarily destabilize the orbit of the body if 
the influence is weak. 

5. Migration to the Main Belt 

About 30 - 40 % of test bodies migrate to the Main Belt on 
non-cometary (asteroidal) orbits (see Hilda type orbit in Fig. 2). 
Most of them leave this region in a dynamically short time, usu- 
ally in less than 250 orbital periods without entering a MMR, 
but about 30 % of them will have resonant interactions with 
Jupiter and/or Saturn. When they enter a MMR, their lifetime 
in the Main Belt grows rapidly - up to 20,000 orbital revolu- 
tions in our calculations. 

The modelling shows that proto-cometary objects in their 
pre-active state can enrich the Main Belt population. A partic- 
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Fig. 2. 1998 SG35 - test particle no 35. A. Variations of the 
semi-major axis as a function of time. The orbit is locked in a 
2/3 MMR with Jupiter after 118,000 years of dynamical evo- 
lution. B. Variations of the critical angle as a function of time. 
When the particle is in resonance, the critical angle starts to li- 
brate about 0° with a large amplitude (dotted-free area indicates 
values of angles which are unaccessible for the librating body). 
The amplitude increases in time and the libration is hardly vis- 
ible at the end of the integration (dotted-free area almost disap- 
pears). The particle will probably leave the resonance in several 
thousand years. When the body remains outside the resonance, 
values of the critical angle vary at random. C. The phase space 
of the eccentricity variations in the function of the semi-major 
axis. When a body is located in a resonance, the eccentricity 
can reach high values but the particle does not leave the MMR. 
The body seems to "stick" to the resonance. 



40,000 80,000 120,000 160,000 200,000 

Time [yr] 

Fig.3. 1999 UGs - test particle no 1. A. During its dynamical 
evolution, the particle enters different MMRs. When in reso- 
nance, the variations of eccentricity (B.) and inclination (C.) 
values evolve slowly (but not periodically) in small boundaries 
due to the lack of close approaches to planets. Most of these 
orbits are Toro class, which means that they cross planetary or- 
bits (in ecliptic projections) but avoid close approaches to these 
massive bodies. 



ularly interesting case is that of the short-lived resonant groups 
of asteroids, like the Griquas. This is a group of resonant aster- 
oids dynamically differentiated from the other resonant groups, 
and that cannot be primordial resonant bodies due to their short 
lifetimes. These vary from 10^ to 10^ years (Roig et al. .2002j 



R. Gabryszewski and I. Wlodarczyk: The Resonant Dynamical Evolution 



5 



Z) 

< 



Close encounters to 
Saturn 




Test particle close to 1/1- and 5/2 MMRs 
with Saturn and Jupiter respectively 



i^P^Nl ill 





Time [yr] 



Fig. 4. 2001 BL41 - test particle no 11 . A. The distance between 
the body and planets as a function of time. B. Variations of 
the semi-major axis as a function of time. The orbit is located 
in overlapping MMRs due to near resonant relations between 
Jupiter and Saturn (so called Great Inequality). Time variations 
of the MMR critical angles of Saturn (C) and Jupiter (D) are 
also presented. After 100,000 years of dynamical evolution, 
when the body remains in 1/1 MMR with Saturn and in 5/2 dis- 
turbing MMR with Jupiter, it has no close approaches to both 
planets. 



which may indicate that a certain per cent of these asteroids 
came from another regions. Have they come from the outer 
Solar System ? Unfortunately, our integrations do not give an 
answer to this question. The calculations indicated only 1 par- 
ticle with 3.1 < a < 3.27 and e > 0.35 close to a 1/2 MMR 
with Jupiter with the orbital lifetime longer than 10'* years. 

We decided to integrate backwards the equations of motion 
of 147 real asteroids moving close to 1/2 MMR with Jupiter 
to find out if these objects are able to come from the region 
outside the orbit of Jupiter. The results shows that 1/2 MMR 
members can be enriched by Centaurs coming from two differ- 
ent regions but the backward simulations cannot be treated as 
the evidence. 



6. Islands of temporary stability 

The numerical experiment allowed one to find temporary "is- 
lands of stability" in the chaotic region of the giant planets - 
bands where particles are trapped for the time estimated for 
lO'* - 10^ years. These "islands" are observed near 6, 6.8, 8.2, 
9.5, 10, 10.5, 11, 11.6, 12.1, 13, 15.1 and 17.3 AU in the giant 
planets region, and also close to 2.8 - 3.1 and 3.4 - 4.2 AU in 
the Main Belt. 

Fig. 5 shows the bands of stability in the {a,e) plane. All 
orbits located in these zones were temporarily stabilized by 
MMRs in the time span of over 10"* years. This period of time 
is comparable with the mean lifetime of the non-resonant par- 
ticles in the region among giant planets estimated by Franklin 
et al. ( 1989 1. Most of the resonances are observed in the Main 
Belt and outside Saturn's orbit. The region between Jupiter and 
Saturn is weakly inhabited. This can be caused by strong gravi- 
tational perturbations from these two planets. The other hardly 
inhabited region (in fact nearly empty) is located outside the 
orbit of Uranus. This can be due to the ensemble selection. 

The most effective MMRs stabilizing dynamics of particles 
are gathered in Table 3. 



7. Discussion and Conclusions 

The region among the giant planets is one of the most chaotic 
in the Solar System. Strong perturbations cause small bodies 
to evolve in the region among giant planets over a dynamically 
short time equal to the 10^ - 10"* years. Our results show that 
this time can be much longer if the body is temporarily trapped 
by a MMR with a planet. 

Temporary stable orbits of the Centaurs seem to be a good 
example of the "resonance stickness" (Dvorak et al. |1998| 
Tsiganis et al. I2000> . Gravitational perturbations can easily 
push a body into a MMR even if the eccentricity of the body 
is high. The evolution of the semi-major axis, eccentricity and 
inclination is chaotic in the borders of the resonance. The per- 
turbations can push the body outside the resonance in timescale 
of 10'* - 10^ years in the same way as they pushed the object 
into this MMR. The number of small bodies which enter a res- 
onance on highly eccentric orbits during a short time of dy- 
namical evolution (200,000 years) indicates that a resonance 
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Table 3. Frequently observed MMRs and their stabilization 
time 



Object 


MMR(s) 


time of orbital 
stabilization 


1998 SG35 


2/3 (J) 


80 kyrs 




2/1 (S) and 5/1 (J) 


50 kyrs 




1/1 (S) and 5/2 (J) 


50 kyrs 




6/7 (U), 1/3 (U), 4/9 (N) 


20 kyrs 


1999 UGi 


2/1 (S) and 5/1 (J) 


60 kyrs 




1/4 (S) and 5/8 (J) 


35 kyrs 




5/4 (S) and 2/9 (N) 


25 kyrs 




1/1 (S) and 5/2 (J) 


20 kyrs 




4/7 (J) 


10 kyrs 


2000 ECgi 


1/4 (S) and 5/8 (J) 


120 kyrs 




5/2 (J), 1/1 (S) and 1/3 (U) 


110 kyrs 




4/9 (J) 


65 kyrs 




1/1 (S) and 5/2 (J) 


30 kyrs 


2001 BUi 


2/1 (S) and 5/1 (J) 


150 kyrs 




2/1 (S), 5/1 (J) and 5/7 (U) 


90 kyrs 




1/1 (S) and 5/2 (J) 


70 kyrs 



sticking is a common mechanism affecting particle lifetimes in 
chaotic zones. 

To ensure that the phenomenon we encountered is in fact 
the resonance stickness, we investigated the autocorrelation 
function of the time series of the osculating elements of one 
of the resonant orbits. If the body is "stuck" to a resonance 
we observed quasi-periodic variations of orbital elements on 
a timescale much longer than the Lyapunov time of this orbit 
(Tsiganis et al. 2000) . 

Fig. 6 presents the autocorrelation function of the semi- 
major axis as a function of the time-lag k for one of the test 
particles being in 1/1 and 5/2 MMRs with Saturn and Jupiter, 
respectively. Time series consisted of 1302 points (which corre- 
sponds to 40,000 years) starting at 1=100,000 years. The series 
is strongly suppressed in time but this is due to the presence of 
gravitational perturbations from the giant planets in our 6 body 
problem. If we had a 3-body problem the suppression would be 
absent. Although the suppression is high, quasi-periodic varia- 
tions can be observed. This seems to confirm the thesis of res- 
onance stickness. 

The most unintuitive result is that overlapping resonances 
can stabilize small body orbits on a much longer timescale than 
single resonances. The overlapping resonances are usually re- 
sponsible for chaotic evolution, but our results show that they 
can stabilize dynamics rather than destabilize it. When we ob- 
serve the particle motion in overlapping resonance with plan- 
ets, we can see that its evolution is chaotic but only inside the 
borders of the MMRs. The resonances protect the body from 
close approaches to planets. Giant planets are close to reso- 
nant relations (Jupiter - Saturn - Uranus and Saturn - Uranus 
- Neptune, Sussmann and Wisdom ll988l[T992t which means 
that a small body can be locked in a MMR with 1 planet and 
close to MMRs with the other 1 or 2 massive bodies. When 
gravitational perturbations change the orbit of the particle, the 
body can be easily "pumped up" into another MMR. Phase 
space of MMRs is thicker when the eccentricity increases and 
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The autocorrelation tunction r(k) ot the semi-major axis as a function of the time-iag k 
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Fig. 6. 2001 BL41 - test particle no 11. The autocorrelation 
function r(k) of the semi-major axis as a function of the time- 
lag k. Time series consists of 1302 points which correspond to 
40,000 years of time. 



hence bodies can remain in overlapping MMRs. The overlap- 
ping resonances protect the body from close approaches to 
planets (which are the main mechanism of chaos in the region) 
and according to our calculations bodies are unable to change 
the portion of the phase space where they evolve, on a timescale 
of over 10^ years. 

Resonances between the test particle and the 4 planets have 
not been observed. This is in good agreement with Sussman's 
and Wisdom's works cited above. 

The results presented in this work indicate that Centaurs 
may also enrich the population of asteroids in the Main Belt. 
MMRs can stabilize Centaur orbits on the timescale of 10"* - 
10^ years. This is a very short period of time in comparison 
with the lifetime of some asteroidal families in this region but 
shows that the Main Belt can include non-processed inactive 
cometary nuclei, not only extinct material. It is hard to say if 
such nuclei can be called primordial but our results show that 
the material can reach the asteroid belt on non-cometary or- 
bits, without approaching the Sun. The particles did not enrich 
specific groups of asteroids in the Main Belt. The only group 
which seems to be overpopulated are Trojans. 

We tried to consider if Centaurs can be a parental popu- 
lation of groups of asteroids with the lifetime comparable to 
10"* - 10^ years. Griquas are an example of such a group. They 
do not originate in the break-up event, a lot of the group mem- 
bers have lifetimes comparable to the lifetimes of Centaurs in 
the Main Belt. But, unfortunately, our modelling did not pro- 
duce a group of objects with 3.1 < a < 3.27 and e > 0.35. 
Although, when we integrated backwards all known orbits of 
the Griqua group members, we found that some of them are 
able to come from the regions among giant planets. Further 
investigations are needed to explain the origin of Griqua type 
bodies. 

The earlier works of Everhart J1973> and Franklin et al. 
(1989) indicate that there are two narrow bands at 7.02 and 
7.54 AU where particles can survive longer than 10^ years. In 
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1997 Holman (I1997I I pointed to a region between 24 and 27 
AU where low inclined (i < 1°) and low eccentric (e < 0.01) 
objects can survive 1 Gyr Our results do not indicate any sta- 
ble orbits close to those bands. This means that the resonant 
dynamics is not the main mechanism responsible for stabiliz- 
ing orbits in these areas. Slow evolution of low inclined, quasi 
circular orbits close to bands at 7.02 and 7.54 AU and in the 
Holman's region is caused by tiny variations of the gravita- 
tional force acting on bodies. 

Centaurs are perceived as a intermediate stage between 
Kuiper Belt objects and Jupiter Family comets (Levison & 
Duncan 1997 1. During the integration we noticed only 10 % 
of orbits with ^ < 1 .5 AU and g > 5.0 AU. Most of them were 
orbits of 1998 5G35 and 2001 BL41 objects (79 and 13 orbits 
respectively). 

It is hard to compare our results with those presented in 
Grazier's and Robutel's papers due to different timescales. 
Their work focused on research of long time stability of Solar 
System regions while we studied the paths of a small ensem- 
ble of orbits. We believe that the results of our experiment and 
Robutel's survey would be similar for the same timescales of 
integrations. 

Further studies are required to broaden our knowledge on 
orbital stabilization by overlapping and many-body MMRs in 
the chaotic region of the Solar System. 
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